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In this paper, we investigate the topological Hubbard model on honeycomb lattice. By considering 
the topological properties of the magnetic state, new types of quantum states - A-type and B-type 
topological spin-density-waves (A-TSDW and B-TSDW) are explored. The low energy physics is 
solely determined by its Chern-Simons-Hopf gauge field theories with different /C-matrices. In the 
formulism of topological field theory, we found spin-charge separated charge-flux binding effect for 
A-TSDW and spin-charge synchronized charge-flux binding effect for B-TSDW. In addition, we 
studied the edge states and quantized Hall effect in different TSDWs. 
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Landau's symmetry breaking paradigm has been very 
successful as a basis for understanding the physics of con- 
ventional solids including metals and (band) insulators. 
In Landau's theory different orders are classified by sym- 
metries. The phase transitions between one type of or- 
dered phase and another one (ordered or disordered) are 
always accompanied by symmetry breaking. Taking spin- 
density-wave (SDW) as an example. To describe such 
ordered state with spontaneous spin rotation symmetry 
breaking, one can define a local order parameter that 
differs for different SDW states, (antiferromagnctic (AF) 
order, ferromagnetic order, ...). 

As the first example beyond the Landau's symmetry 
breaking paradigm, the integer quantum Hall (IQH) ef- 
fect is a remarkable achievement in condensed matter 
physics[l|. To describe the IQH state, the Chcrn number 
or so called TKNN number, Q, is introduced by inte- 
grating over the Brillouin zone (BZ) of the Berry field 
strength [2| . Another example is the fractional quantum 
Hall (FQH) effect 0,13, whose elementary excitations are 
anyons and the subtle structure is dubbed the topological 
order @, HI. Recently, a new class of topological state 
- topological insulator (TI) is found with the quantized 
spin Hall effect 0, Q. For all these topological states, ele- 
mentary excitations are gapped. Particularly, there is no 
local order parameter to characterize them. Instead, the 
low energy properties of these topological states can be 
described by effective Chern-Simons (CS) theories. 

In this paper, we will focus on a new type of topological 
quantum states with spontaneous spin rotation symme- 
try breaking - topological SDW (TSDW) states. Differ- 
ent TSDW states have identical local order parameter - 
the staggered magnetization. In order to further classify 
them, we derive effective Chern-Simons-Hopf (CSH) the- 
ories with different topological matrices, /C-matriccs that 
have been introduced in FQH states @. 

The topological Hubbard model: The Hamiltonian of 
the topological Hubbard model on honeycomb lattice is 
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Here Hh is the Hamiltonian of Haldane model (l2j 
which is given by Hn = -t ^ {c\ a c J(J + h.c.) - 

t' e^^c^Cjcr. t and i' are the nearest neighbor and 

the next nearest neighbor hoppings, respectively. We in- 
troduce a complex phase 0y (\<fiij \ = to the next near- 
est neighbor hopping, of which the positive phase is set to 
be clockwise. H ' = e clavier — £ ^Lr^Ur denotes 

an on-site staggered energy : e on A site and — e on B 
site. U is the on-site Coulomb repulsion, fj. is the chemi- 
cal potential and fj. = U/2 at half-filling for our concern 
in this paper. 

When U is zero, the spectrum for free 
fermions is = ±yj |£fc| 2 + + e) 2 where 

|£ k | = + 2 cos (VSk y ) + 4 cos (3^/2) cos (V3k y /2) 

and C k = 2t ' [—2 cos (3*^/2) sin(v / 3fc w /2) + sin(-\/3fej,)] . 
According to this spectrum, we can see that there exist 
energy gaps A/i, A/ 2 near points ki = — ^(1, l/v3) 
and k 2 = ^(1, 1/V3) as A fl = \2s-6y/3t'\ and 
A/2 = 2e + 6y/3t', respectively. In addition, there exist 
two phases for this case, the quantum anomalous Hall 
(QAH) state and the normal band insulator (NI) state 
with trivial topological properties. They are separated 
by the phase boundary A/i =0. In the QAH state, due 
to the nonzero TKNN number, Q = 2, there exists the 
IQH effect with a quantized (charge) Hall conductivity 
(T H = 2e 2 /h. 

Mean field (MF) phase diagram and topological quan- 
tum phase transitions: Turning on the interaction, the 
topological Hubbard model is unstable against AF SDW 
order that is described by {c\ a Ci^) = j(\ + (— 1)VM) 
where the local order parameter M is the staggered mag- 
netization. We set a = +1 for spin-up and a = — 1 for 
spin-down, and then for M^O, the Hamiltonian can be 
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written as H = i?n + H' — (~^T Aa/c^ct^c^, with 

Am = UM/2. By MF approach, we obtain the self- 
consistency equation for M by minimizing the energy at 
zero temperature in the reduced BZ as 
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where N s 



2E kl 
the number 



2Eu 



is the number of unit 
v/^+Am+e^ + I&I 2 and E k 



(2) 



cells, 



y/(? k -A M +e)* + \€ k \2. 

To determine the phase diagram, there are two types 
of phase transitions : one is the magnetic phase transi- 
tion [denoted by ([//^)m] between a magnetic order state 
with M ^ and a non-magnetic state with M = 0, 
the other one is the topological phase transition that 
is characterized by the condition of zero fermion's en- 
ergy gaps, A/i = -6a/3£' + 2e + 2A M = or A /2 = 
6\/3t' + 2e — 2Am = 0. After determining the phase 
boundaries, we plot the phase diagram with five different 
quantum phases in FIG.l for e = 0.15: QAH state, NI, 
A-typc topological SDW state (A-TSDW), B-typc topo- 
logical SDW state (B-TSDW), and trivial SDW state. In 
FIG. 2, we also plot the staggered magnetization and the 
energy gaps A/i, A/2 for the same case. 
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FIG. 1: (Color online) The phase diagram for the case of 
e = 0.15 : I is QAH state, II is A-TSDW, III is B-TSDW, IV 
is trivial SDW, V is NI. The black, red and blue lines are the 
critical lines of (U/t)M, (U/t) c i and (U/t) C 2, respectively. 

Based on the MF results, the TKNN numbers in A- 
TSDW, B-TSDW and the trivial SDW states are Q = 2, 
Q = 1, Q = 0, respectively. However, due to the quan- 
tum spin fluctuations, the classification of SDW states by 
the TKNN number is insufficient to give the final answer. 
Instead, a 2-by-2 matrix (/C-matrix) plays a key role in 
the topological classification of the SDW orders with the 
same local order parameter, M. 

Induced CSH terms and K -matrices representation of 
TSDWs: In this part we will derive the low energy effec- 
tive theory of (T-)SDW states by considering quantum 
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FIG. 2: (Color online) The staggered magnetization and the 
energy gaps A/i, A/2 for the fermion excitations at node ki, 
k2 for the case of e = 0.15 and t' = 0.1. Below {U/i)M, 
A/i = A/2, that is the black line; above (U/t)u, the red line 
denotes A/i and the blue line denotes A/2. 



fluctuations of effective spin moments based on a for- 
mulation by keeping spin rotation symmetry, o~ z — > n • a 

where n is the SDW order parameter, (c\aci) = Mn. In 

this case, the Dirac-like effective Lagrangian with spin ro- 
tation symmetry in the continuum limit can be obtained 
as 

Cf = [iipa'Yp (<9 M - iAp) ip a + m a ^ a 4' a - SA M ip a a ■ nip a ] 

a 

. (3) 

which describes low energy charged fermionic 
modes a = 1 near ki, 4>i = V4to = 
( ifoiA, i>fiB, 4>iia, i'liB ) and a = 2 near k 2 , 

V>2 = IptlO = ( i'f2B, $f2A, 1pl2B, 1pl2A )■ TllC 

masses of two-flavor fermions are mi = e — 

and m% = e + 3\/3Y. 7^ is defined as 70 = cro ® i~ z , 

7l = CT <8> T y , 72 = O-q <g) T x With (7 = ( q j Y T x , T y , 

t z are Pauli matrices. 5—1 for a = 1 and 5 = — 1 for 
a = 2. We have set the Fermi velocity to be unit Vp = 1. 

In CP 1 representation, we may rewrite the effective 
Lagrangian of fermions in Eq. ([3]) as 

£f = E ^ i^ 9 ^ + 7m^m _ 7^3^ + m a - 6A M o- 3 ) ip' a 

a 

(4) 

with ip a (r, r) = U' (r, r) ip a (r, r) , where U (r, r) is a local 
and time-dependent spin SU(2) transformation defined 
by W (r, r) n • aU (r,r) = 03. And a M is introduced as 
an assistant gauge field as 103 a M = W (r, r) c^f/ (r, t) . 

An important property of above model in Eq.Q is the 
current anomaly. The vacuum expectation value of the 
fermionic current J£ a = i^a^l^^Pa.a) can be defined by 
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Jta = *{7"[(*-D + im a , a )\it) + im^)}' 1 (it) + im a ^} 
where D — 7 M (9 M — iA^ + iaa^) and the mass terms are 
ni>a,cr = nT>a — SA^ja. The topological current J£ a is 
obtained to be 

Then we derive the CSH terms as Ccsh = —i J2a ct(^a« — 

To make an explicit description of TSDWs, we intro- 
duce the /C-matrix formulation that has been used to 
characterize FQH fluids successfully Q. Now the CSH 
term is written as 

£ ^ = -«E^^<3X (6) 
i, j 

where JC is 2-by-2 matrix, a^ =1 = and a^ =2 = a^- 
The " charge" of A^ and a M are defined by q and q s , re- 
spectively. 

Thus for different SDW orders with the same order pa- 
rameter M, we have different /C-matrices : for A-TSDW 

order with mi, m 2 > A M , JC = \ l ° J ; for B -TSDW 

order with ni2 > Am > mi, JC = ( j 1 J ' ^ v * a ^ 

SDW order with mi, m2 < A a/, /C = 0. It is obvious that 
such topological structure labeled by /C-matrices is be- 
yond Landau's symmetry breaking paradigm and TKNN 
number classification. 

In the following parts we will use the following effec- 
tive model with the CSH term to learn the topological 
properties of different SDW orders, £ c ff = Cf + Ccsh- 

A-TSDW: In A-TSDW, an important property is 
" spin- charge separated charge-flux binding" effect for 
gauge fields a\ and A\. From the effective CSH La- 
grangian in Eq.©, we get the equations of motion for 
a\ and A\, 

-e^duAx = - J M , -e^ x d v a x = -J s>1 (7) 
where _J M = i(J2 a ^al^a) and J sll = 

-iiYta^'al^^'a) = ~*(Ea ^al^ ■ CT^a) ■ As a re- 
Suit, we get the identities $ e /7r = — q and § s /tt = —q s . 

On the one hand, the unit magnetic flux of A\, $ e = 
J d 2 re°" x dvA\, will bind electric charge q = —2. For ex- 
ample, a half magnetic flux $ e = ir will carry q = — 1 
electric charge at its core (see FIG. 3. (a)). On the other 
hand, the unit "magnetic flux" of a\ will bind "charge" 
q s = -2 (sec FIG.3.(b)). The "magnetic flux" of a M 
is in fact the topological spin texture of SDW order (so 
called Q skyrmion) that is characterized by winding num- 
ber Q = j- J d 2 reoi/An • 9„n x (9\n[ll[ from the relation 
d u a\ — d\a v = \n ■ d u n x d\n. Thus in A-TSDW, each 
skyrmion with <£> s = 2irQ "magnetic flux" will accom- 
pany by additional "charge" number q s = —2Q. As a 
result, skyrmion's spin is S = Q 2 <d/2n with 6 = 2n. For 
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FIG. 3: (Color online) Spin-charge separated charge-flux 
binding effect of A-TSDW : the induced quantum numbers 
on a half magnetic flux (a) and those on a Q = 1 skyrmion 
(b). 



example, by binding q s = —2 "charge" , the spin of Q = 1 
skyrmion is S = I. 

Furthermore, if there exists a tiny easy-plane 
anisotropic term (for example, the spin-orbital coupling 
term), we can define the skyrmion with fractional wind- 
ing number - the half skyrmion with Q = ±1/2, of which 
the induced "charge" is g s = =pl. However, because the 
energy of the half skyrmion in a long range SDW or- 
der diverges logarithmically, we cannot treat it as a real 
quasiparticle. 

In A-TSDW, from the CS term - ^e^ x A^d u A x and 
Hopf term —-^c^^a^d^ax, we find four right-moving 
branches of edge excitations instead of two, which are 
described by the following one dimension (ID) fcrmion 
theory 

a 

(8) 

where a, j3 = 1,2. ip c ^ a carries a unit of A^ charge and 
ips,p a urr it °f a [i chargeQ. That means we get spin- 
charge separated edge states : two edge modes only carry 
electric charge current; two only carry spin current. 

Consequently, we can define two types of Hall 
conductivities : the quantized charge Hall con- 
ductivity <j H = lim^^o (Ji(uJ,0)J j (-uj,0)) 
and the quantized spin Hall conductivity er s = 
lim w _>.o ~ey {J si (u;,0)J s j(— w,0)) . Here J ?; denotes 
electric current, Ji = i(%2 a V'aTiV'a); and J S i denotes 
spin current, J s i = —i(J2 a V'aTi 11 ' o"ipa)- From the 
CSH term, we obtain quantized Hall conductivities : 
<jh = a 3 = 2e 2 /h that correspond to the spin-charge 
separated edge states [Toj. 

B-TSDW: In B-TSDW, the equations of motion for a p 
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and turn into 

±e» vX d v A x + -U^ A cW = - J M (9) 

Z7T Z7T 

and 

J- e ^^^ A + -U^ A cW = - J SM . (10) 

Z7T Z7T 

As a result, in B-TSDW, due to " spin- charge syn- 
chronized charge-flux binding" effect from the identities 
<J> s /27r + <I> e /27r = — q = —q 8 , the induced electric charge 
number is always equal to the induced "charge" number 
of a M on a topological object. 
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FIG. 4: (Color online) Spin-charge synchronized charge-flux 
binding effect of B-TSDW : the induced quantum numbers on 
a Q — 1 skyrmion (a) and those on a half magnetic flux (b). 

On the one hand, due to the condition $ s /27r = — q = 
— q s , Q skyrmion will carry q = —Q electric charge of 
gauge field and q s = —Q "charge" of gauge field 
a M . For example, Q = 1 skyrmion carries a unit elec- 
tric charge q = — 1 and a unit "charge" q s = — 1 (see 
FIG.4.(a)). With a unit "charge" q s , Q = 1 skyrmion 
gets half spin and becomes a charged 5=1/2 fermion; 
For a Q — 2 skyrmion, there exist two induced charge 
numbers on it, q = —2, q s = —2. Then it becomes 
a charged 5 = 2 boson. In addition, for a SDW or- 
der with easy-plane anisotropic energy, we have the half 
skyrmion with fractional winding number Q = ±1/2, of 
which there exist fractional charge numbers, q = 



q s = Tl/2. On the other hand, using the same approach, 
from $ e /27r = — q = —q s , we find that a half quantized 
magnetic flux $ e = tt will also carry fractional charge 
numbers, q = -1/2, q s = -1/2 (see FIG.4.(b)). 

For B-TSDW, due to a mutual CS term, 
— j^e^ x Afj,d v a\, there is no spin-charge separated 
edge states. Instead, we have spin-charge synchronized 
edge states. Now the charge of a+ lfl = A^ + 
is quantized to an integer number. Then the ef- 
fective CSH field theory has one right-moving edge 
excitation. The edge excitation is described by the 
following ID fermion theory £ d g c = ip'(dt — vd x )ip 
where -0 carries a unit of a+, M " charge" 0. Then wc 
may define the spin-charge synchronized Hall conduc- 
tivity as o = lim w _^o zi e ij ^i( w : 0) Jj(— w, 0)^ with 

Ji = i(E a fa,^fa,i) = *(Ea ^-7<(l ~ n • *)i> a /2) 

and get a quantized spin-charge synchronized Hall 

2 

conductivity : a = 2_ that corresponds to the edge 
state. That means an electric field can drive both a 
quantized electric charge current and a quantized spin 
current. In particular, such a quantized spin-charge 
synchronized Hall effect is not QAH effect for electrons 
with confined spin-charge degrees of freedom in QAH 
state. 

Conclusion: In this paper, we investigate the topolog- 
ical Hubbard model on honeycomb lattice. New types 
of quantum states - A-TSDW and B-TSDW states are 
explored which bear an identical staggered magnetiza- 
tion M as the local order parameter. To characterize 
different TSDWs, we introduce ^C-matrices to denote dif- 
ferent CSH gauge theories. In the formulism of the CSH 
gauge field theories with different ^-matrices, we found 
spin-charge separated charge-flux binding effect in A- 
TSDW and spin-charge synchronized charge-flux binding 
effect in B-TSDW. In addition wc studies the edge states 
and corresponding quantized Hall effect in TSDWs. Our 
findings suggest that although there exists spontaneous 
spin rotation symmetry breaking, the TSDWs are beyond 
Landau's paradigm. 
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A. The detailed calculations of the 
Chern-Simons-Hopf (CSH) terms in Eq.(6) 

Firstly we calculate the induced Chcrn-Simons (CS) 
term of a one flavor fermionic-er model. The Lagrangian 
of one flavor fcrmionic-er model is written as 

£ = #7/i(fyt ~ ibfj,)4> + rmpip 

where to is a fcrmion mass. To obtain the induced CS 
term, we integrating over fermions and get 

S e s = hi Z 

where 

Z = J [dip] [dip] exp (- J dxC). 

The one fermion loop effective action becomes 
S e ff = lndet(i7 M 3 A1 + 7 M 6 M + m) 



(11) 



trlog(i7 M <9 M + m) + tr( 
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1 / 1 

7: tr r 

2 «7/ 1 o ft + to. 



7/ifyx" — 7/A) 

n^o^ + to 



Then the quadratic term of in the effective action is 



S cS = - 



d 3 p 
(2tt) 3 



[b^(-p)D^b^p)} (12) 



where D^ u is 
d 3 k 



(2tt)= 



tr 



M £V7/x + V7ji - m v k } a E --rn 
7 /_ , 1 \o 1 9 7 ;o . o 



(p + k) 2 + to 2 

Under trf^^ 1 ^*) = -2e^ A , we obtain 

d 3 fc 1 



(13) 



D» v (p,m) = e^ vX p x 2m 



e^ X Px 



(27r) 3 [(p + k) 2 + TO 2 ] [fc 2 + TO 2 ] 

(14) 

1 TO . , \p\ 



arcsin(- 
2tt \p\ ^/p2 + 4 TO 2 



In the long wavelength limit, (^ — > 0), due to O 
SrTmT' we an m duced CS term as 



8n to 



(15) 



Next we calculate the induced CSH term in Eq.(6) in 
the paper that denotes a four-component fermionic model 

At = V'i-r + 7 /1 ^ ( u - 7% + ™i - a m) V'if 

£14 = #4. + 'fAp + 7% + to! + A M ) ^4- 

-Cat = V^2f (*7 M ^ + 7^4 _ 7^ + to 2 + A M ) 

£24 = 4>2i (h^n + 7^ + 7% + ?7i 2 - Am) ^24- 
Integrating we get the induced CS term 

87T |mi - Am I 
Integrating ip^, , we get the induced CS term 



3tt Itoi + A 



Ml 

Integrating ^-fi wc get the induced CS term 

C n {A») = -± ™ 2 + ^ M e^\A, - a,)d v {A x - 
87r |to 2 + Am I 

Integrating tj)^, we get the induced CS term 

C 2i (A,) = 4 r"A M ^(^ + «^(A A + a A ). 
8tt |to 2 - Am I 

Then we get the total induced CSH term as 

C(A^) = £i t (Ai) + ZuiAv) + £»t(4*) + 
(i) For the case of mi, to 2 > Am, we get 

£(-A„) = £i t (4J + £ u (4) + £ 2t (4) + £ 2 ^(A M ) 



7 

8tt 



— e^ x (A, l + a^(A x + a x ) 



e^ x (A^- afl )d v (A x -a x ) 



-—e^ x (A^-a^(A x -a x ) 



8- 
/ 

Sn 
i 



—e^ x (A fl + a f ,)d l/ (A x + a x ) 



57T 



(ii) For the case of to 2 > A a/ > mi, we get 

£(A M ) = £i t (4) + £u(4) + £2t(4) + £2i(4) 
= l _ f ^x 



87T 



2 

tin 



+ — c^ A (A M -a M )ft,(A A -o A ) 



, e» vX {A„-a»)d v {A x -a x ) 

tin 

^-^(Ap + aJdviAx + ax) 

tin 



-in 



\» vX A il d v Ax + ^ x A,d v a x + -le^ x a,d„a x - 



2tt 



47T 



6 



(iii) For the case of mi , rri2 < Am , we get Finally, by the /C-matrix, the CSH term is written as 

UA r ) = C„(A r ) + C u (A f ) + C 2 ,(A r ) + C u (A r ) 

= '^ (All + « MAx+n) C CSH = -iY,-£ c »^ai ,16) 

■ where JC is 2-by-2 matrix, a^ =1 = A^ and a^ =2 = a^. 

— — e M — a^)d v {A\ — a\) Thus for different SDW orders with the same order pa- 
rameter M, wc have different /C-matrices : for mi, m-2 > 

+ + <*>*^ + ^ A M , K = (I I ) ; for m 2 > Am > m lt K = ( \ \ 



0. 



for mi, m-i < Am, JC = 0. 



